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ABSTRACT 

We present high angular resolution observations (0.5" xO.3") carried out with the Submillimeter Array 
(SMA) toward the AFGL2591 high-mass star forming region. Our SMA images reveal a clear chemical 
segregation within the AFGL2591 VLA 3 hot core, where different molecular species (Type I, II and 
III) appear distributed in three concentric shells. This is the first time that such a chemical segregation 
is ever reported at linear scales <3000AU within a hot core. While Type I species (H2S and ^^CS) 
peak at the AFGL2591 VLA 3 protostar, Type II molecules (HC3N, OCS, SO and SO2) show a double- 
peaked structure circumventing the continuum peak. Type HI species, represented by CH3OH, form 
a ring-like structure surrounding the continuum emission. The excitation temperatures of SO2, HC3N 
and CH3OH (185±11K, 150±20K and 124±12K, respectively) show a temperature gradient within 
the AFGL2591 VLA 3 envelope, consistent with previous observations and modeling of the source. 
By combining the H2S, SO2 and CH3OH images, representative of the three concentric shells, we 
find that the global kinematics of the molecular gas follow Keplerian-like rotation around a 40 M©- 
star. The chemical segregation observed toward AFGL2591 VLA 3 is explained by the combination 
of molecular UV photo-dissociation and a high-temperature (~1000 K) gas-phase chemistry within 
the low extinction innermost region in the AFGL2591 VLA 3 hot core. 

Subject headings: stars: formation — ISM: individual objects (AFGL2591) — ISM: molecules 



1. INTRODUCTION 

Hot cores are hot (~200K), compact (<0.1pc) and 
dense (>10^ cm~^) condensations, which represent one of 
the earliest stages of massive star formation (e.g. Garay 
& Lizano 1999). These objects are chemically very rich, 
and show large enhancements of S-bearing species such 
as H2S, OCS or SO2, and of complex organic molecules 
(COMs) such as H2CO, CH3OH or CH3CN (Blake et al. 
1987; Comito et al. 2005; Hatcheh et al. 1998b). These 
enhancements are attributed to the thermal evaporation 
of the mantles of dust grains, followed by rapid ion- 
neutral gas-phase reactions (Charnley et al. 1992; Charn- 
ley 1997; Hatcheh et al. 1998a; Viti et al. 2004; Wakelam 
et al. 2004). 

Toward hot cores, the warming up of the gas in the 
envelope by the central protostar, generates a tempera- 
ture gradient as a function of radius that largely affects 
the chemical composition of the grain mantles and of the 
molecular gas within these cores. As a consequence, a 
chemical segregation is expected to occur (van Dishoeck 
& Blake 1998; Nomura & Mihar 2004). Indirect evidence 
of this molecular segregation has been found in the form 
of abundance discontinuities for COMs toward hot cores 
(van der Tak et al. 1999, 2000b), as well as toward low- 
mass warm cores (or hot corinos; see Ceccarelli et al. 
2000; Maret et al. 2004; Awad et al. 2010). However, 
direct imaging of the chemical segregation at the inner 
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regions of hot molecular cores, still remains to be re- 
ported. 

The AFGL2591 high-mass star forming region is lo- 
cated in the Cygnus X complex at a distance of ~3 kpc. 
This distance, which was substantially underestimated 
in the past, has been derived recently by using trigono- 
metric parallax of the 22.2 GHz H2O masers detected 
toward this region (Rygl et al. 2012). The revised mass 
of the large-scale AFGL2591 clump is 2xlO'^M0, and 
the measured total IR luminosity is ~2xlO^L0 (Sanna 
et al. 2012). AFGL2591 hosts a cluster of young stel- 
lar objects (YSOs) over scales of 0.1 pc (Campbell 1984; 
Trinidad et al. 2003; van der Tak & Menten 2005), where 
several radio continuum sources are detected (see VLA 
1, VLA 2 and VLA 3 in Trinidad et al. 2003). Sanna 
et al. (2012) has recently reported the detection of an- 
other YSO (the NE source) at -0.4" (1300 AU) north 
the AFGL2591 VLA 3 object. 

From ah the YSOs detected in the AFGL2591 star 
forming region, AFGL2591 VLA 3 is believed to be the 
youngest (dynamical age of -2xl0'^yr; Doty et al. 2002; 
Stauber et al. 2005) and most massive source in the 
cluster (estimated mass of — 38M0; Sanna et al. 2012). 
AFGL2591 VLA 3 is very bright in the mid-IR (de Wit 
et al. 2009), and is responsible for most of the IR lumi- 
nosity measured toward this region (Sanna et al. 2012). 
AFGL2591 VLA 3 is also known to power a large-scale 
east- west outflow detected in CO and near-IR emission 
(P.A.-280°; Mitchell et al. 1992; Preibisch et al. 2003), 
whose axis lies close to the direction of the line-of-sight 
(Hasegawa & Mitcheh 1995; van der Tak et al. 1999). 
The blue-shifted lobe of the outflow is detected toward 
the west of AFGL2591 VLA 3, while the red-shifted out- 
flowing gas is found toward the east (Mitchell et al. 1992). 
The AFGL2591 VLA 3 molecular envelope is seen almost 
pole-on, and rotates in the counter-clockwise direction 
(van der Tak et al. 2006). van der Tak et al. (1999) have 
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proposed that the AFGL2591 VLA 3 molecular envelope 
has two different physical regimes with an inner and hot- 
ter core and an outer and cooler envelope. This makes 
the AFGL2591 VLA 3 hot core an excellent candidate to 
probe for chemical segregation. 

In this paper, we report the first detection of a clear 
chemical segregation at linear scales <3000AU within 
the hot core around the AFGL2591 VLA 3 high-mass 
protostar. The comparison of the SMA images with gas- 
grain chemical models shows that the observed chem- 
ical segregation is due to i) a strong UV radiation 
field leading to molecular photo-dissociation; and ii) a 
high-temperature gas-phase chemistry within the low- 
extinction innermost region in this core. The paper is 
organized as follows. The SMA observations are reported 
in Section 2. The images of the continuum and molecu- 
lar line emission toward AFGL2591 VLA 3 are presented 
in Section 3. The derived excitation temperatures, col- 
umn densities and abundances of the molecular species 
measured toward this source are given in Section 4. And 
the modelling of the chemical segregation detected to- 
ward this object is presented in Section 5. In Section 6, 
we summarize our conclusions. 

2. OBSERVATIONS 

Observations of the AFGL2591 star forming region 
were carried out with the SMA^ in the very extended 
configuration (VEX) on 2010 February 15. The size 
of the synthesized beam of our VEX observations was 
0.55" X 0.34" (see Figure 1 for the UV coverage and 
dirty beam of the VEX data). The phase center of 
the observations was set at a( J2000)=20^29^24.90^ 
(5(J2000)=+40°11'20.3''. The zenith opacity at 225 GHz 
was 0.06 and the double sideband system temperatures 
were typically 200-280 K. The receivers were tuned to 
an LO frequency of 218.75 GHz, and the correlator pro- 
vided a uniform spectral resolution of 0.8 MHz (i.e. 
~l.lkms~^). The radial velocity of the source was set 
at -5.5kms-i (van der Tak et al. 1999). We used 3C273 
(12.5 Jy) as bandpass calibrator, and MWC349A (1.7 Jy) 
as flux and gain calibrator. Calibration of the raw data 
was carried out within the IDL MIR software package, 
and continuum subtraction, imaging and deconvolution 
was done within MIRIAD. 

In order to quantitatively measure the effect of miss- 
ing fiux in our VEX data (see Section 3.4), we have 
also used unpublished data obtained by us with the 
SMA in Subcompact configuration (beam of 6.7" x5.6"). 
The central coordinates and spectral resolution of the 
SMA Subcompact images are the same as those of the 
VEX observations. However, the frequency coverage was 
slightly different from that of the VEX data, and the H2S 
(22,o-^2i^i) hue was not observed. A detailed analysis of 
the large-scale structure of the continuum and molecu- 
lar line emission toward AFGL2591 VLA 3 will be pre- 
sented elsewhere (Jimenez-Serra et al. in preparation). 
In the present paper, we only use the Subcompact data 
to demonstrate that the chemical segregation detected 
toward this source is not an effect of missing short spac- 

^ The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics, and is funded by the Smith- 
sonian Institution and the Academia Sinica. 
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Fig. 1. — UV coverage (upper panel) and dirty beam (lower 
panel) of our observations taken with the SMA in VEX configura- 
tion. The synthesized beam is 0.55" x 0.34" (Section 2). 



ings in the VEX observations. 

3. RESULTS 

3.1. Dust continuum emission and ^^CO outflowing gas 

In Figure 2, we show the ^^CO images for the red- 
and blue-shifted gas detected toward AFGL2591 VLA 
3 (red and blue contours, respectively), superimposed 
on the 1.3 mm dust continuum emission (gray scale 
and black contours). The continuum emission (size of 
~0.6"x0.45") is associated with the brightest mid-IR 
source in the region (de Wit et al. 2009), and peaks 
at a(J2000)=20^29^24.889^ ^(J2000)=+40°ln9.5r^ 
This position is very similar to that reported by Benz 
et al. (2007) for the continuum source seen at 348 GHz 
in their SMA images. The derived integrated 1.3 mm 
continuum flux (^lOOmJy) is smaller than previously 
reported ('^190mJy; van der Tak et al. 2006), because 
the PdBI observations (beam of ~0.9") were more sen- 
sitive to extended emission. The VLA 3 source detected 
at 3.6 cm (Trinidad et al. 2003) is located at ~0.2" east 
of the 1.3 mm dust continuum peak, possibly linked to 
the blue-shifted part of the high-velocity jet (Sanna et 
al. 2012). The fainter 1.3mm continuum component to- 
ward the south of AFGL2591 VLA 3 could be associated 
with shocked gas at the walls of the outflow cavity (see 
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Fig. 2. — Integrated intensity maps of the ■'^^CO J=2— )►! emission 
detected toward AFGL2591 VLA 3 from -39.5 to -23kms-i (blue 
contours) and from 8.9 to 38.6 kms"-*^ (red contours), superimposed 
on the 1.3 mm dust continuum map (gray scale and black contours). 
The contour levels for the 1.3mm emission are 2.4 (3cr), 4.0, 5.6, 
17.6, 29.6, 41.6 and 53.6 mJy/beam. For ^^CO, the first contour 
and step levels are 0.30 (2cr) and 0.15 Jy/beam km s"-*^ for the blue- 
shifted lobe, and 0.36 (2cr) and 0.54 Jy /beam kms "-"^ for the red- 
shifted lobe. Negative contours correspond to the 3cr level. Filled 
star and filled square show the location of the VLA 3 and NE 
sources (Trinidad et al. 2003; Sanna et al. 2012). The beam is 
shown at the lower right corner. 



Maury et al. 2010). 

Toward the north of VLA 3 (see Figure 2), the 1.3 mm 
continuum emission shows an elongation that could be 
associated with the NE source reported by Sanna et al. 
(2012). This source, which is not detected at centimeter 
wavelengths (Trinidad et al. 2003), has likely generated 
the double bow-shock structure seen in H2O masers at 
0.4" north VLA 3. As shown in Sections 3.2 and 3.3, the 
NE source does not have any clear counterpart in molec- 
ular line emission. The kinematics of the molecular gas 
toward AFGL2591 VLA 3 seem to be dominated by the 
kinematics of the VLA 3 molecular envelope, although 
the angular resolution of our observations may not be 
sufficient to resolve any molecular feature arising from 
the NE source. 

From Figure 2, we find that the ^^CO emission forms a 
bi-conical structure aligned in the east- west direction and 
centered at the 1.3 mm VLA 3 source. The kinematics 
of this structure are consistent with those of the large- 
scale ^^CO outflow (blue-shifted gas toward the west and 
red-shifted emission toward the east; see Mitchell et al. 
1992). In addition, the ^^CO blue lobe detected in our 
VEX images matches the morphology of the bright near- 
IR loop detected in the K-band (Preibisch et al. 2003). 
The opening angle of the bi-conical structure detected 
with the SMA is --90°-100°, similar to that measured 
from H2O masers (110°; Sanna et al. 2012) or from near- 
IR emission (>100°; Preibisch et al. 2003). 
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3.2. Differences in the morphology of the molecular line 
emission toward AFGL2591 VLA 3 

In Table 1, we report the molecular rotational lines 
measured within the 8 GHz passband of the SMA, and 
Figure 3 presents the integrated intensity maps for some 
of these lines superimposed on the 1.3 mm continuum 
image of the AFGL2591 VLA 3 source. From Figure 3, 
we find that the observed molecular species show differ- 
ent spatial distributions in a clear chemical segregation. 
According to the observed morphology, we classify these 
molecules into Type I, II and III species (see Figure 3). 
While Type I species (H2S) have a compact structure 
peaking at the AFGL2591 VLA 3 continuum emission 
(left panels). Type II molecules (HC3N, OCS and SO2) 
show a double-peaked morphology that circumvents the 
continuum peak (middle panels). For Type III species 
(represented by CII3OH; right panels), the emission is 
distributed in a ring-like structure surrounding the con- 
tinuum core. 

Two species, namely ^^CS and SO, have been ten- 
tatively classified as Type I and Type II species, re- 
spectively. For ^^CS, although the brightest emission 
is detected toward the continuum peak, this molecule 
also presents a secondary peak toward the south-west of 
AFGL2591 VLA 3. This peak does not coincide with any 
of the sources detected in the region or with the molec- 
ular emission peaks from Type II or III species. The 
secondary peak could be associated with outfiowing gas. 
However, its kinematics (at red-shifted velocities) are not 
consistent with those of the outfiow toward this position. 
As shown in Section 5, CS is also expected to be relatively 
abundant toward the regions where Type II species are 
detected, suggesting that the CS emission could be more 
extended than detected with the rarer ^^CS isotopologue. 
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Fig. 3. — Integrated intensity images of several molecular species from -8.7 to -2.1 kms"-*^ (thick contours), superimposed on the 1.3mm 
continuum map (color scale and dashed contour). Filled star, filled square and filled circle show, respectively, the location of the VLA 
3 source (Trinidad et al. 2003), the NE source (Sanna et al. 2012), and the main 1.3 mm dust continuum peak detected with the SMA 
(Section 3.1). The molecular species, transition and energy of the upper level of every map are shown in every panel. The first contour and 
step levels are: 0.27 (3cr) and 0.45 Jy/beam km s-^ for H2S (22,o^2i,i), 0.24 (3cr) and 0.24 Jy/beam km s-^ for ^^CS (5^4), 0.21 (3cr) and 
0.35 Jy/beamkms-i for HC3N i;=0 (24^23), 0.21 {3a) and 0.21Jy/beam km s"! for OCS (18^17), 0.21 {3a) and 0.28 Jy/beam km s"! for 
SO2 (115,7^124,8), 0.36 (3c7) andO.60 Jy/beamkms-i for SO (56^45), 0.21 {3a) and 0.21Jy/beam km s-^ for CH3OH vt=0 (4+2,2^3+1,2 
E), 0.18 {2a) and 0.18 Jy/beam km s-^ for CH3OH vt=0 (102,8^93,7 A+), 0.15 {2a) and 0.15 Jy/beam km s-^ for CH3OH vt=l (6i,5^72'6 
A~), and 0.14 (2c7) and 0.07 Jy/beam km for CH3OH vt=0 (20+i, 19^20+0,20 E). Negative contours correspond to the 3c7 level. Dashed 
lines show the direction of the PV cuts reported in Figure 5. Tl and T2 are the brightest emission peaks seen in SO2 and CH3OH, where 
Tex has been estimated (Section 4 and Figure 7). Beams are shown at the lower right corner in every panel. 



For SO, the double-peaked morphology in the SO inte- 
grated intensity map of Figure 3, is not as clearly seen as 
for HC3N, OCS or SO2. However, a closer inspection of 
the Position- Velocity diagram of the SO emission reveals 
a similar kinematic structure to that of e.g. SO2 (see 
Figures and below). Nevertheless, we cannot rule out 
the possibility that the emission from the low-excitation 
SO (56^45) transition is affected by large optical depth 
effects and/or it is part of the base of the AFGL2591 
VLA 3 outflow. A significant amount of SO could also 
be present within the same region where Type I species 
arise (i.e. H2S and ^^CS), as suggested by the model- 
ing of the chemistry of the AFGL2591 VLA 3 source 
(Sections). 

We note that the double-peaked structure of Type II 
species such as SO2 or OCS, and the ring- like morphol- 
ogy of CH3OH, are not likely produced by excitation 
or optical depth effects since higher-excitation, optically 
thin transitions such as SO2 (227^x5^236, is) (£^n=353K) 
or CH3OH vt=l (61,5^72,6 A-) and CH3OH ^^=0 
(20+1,19^20+0,20 E) (K=374 K and 509 K, respectively) 
present similar spatial distribution and kinematic struc- 
ture to those of the lower-excitation lines (see Figure 3, 
Figures and Section 3.3 below). This behavior clearly 
contrasts with that observed toward the G29.96-0.02 
hot core (Beuther et al. 2007), where the CH3OH high- 
excitation, optically thin lines peak toward the hot core 
while the low-excitation, optically thick transitions show 
a double-peaked morphology surrounding the core, as 
expected from large optical depths and if the CH3OH 
gas were evenly distributed across the hot core. Al- 



though chemical differentiation has been reported toward 
G29.96-0.02 (Beuther et al. 2007), this differentiation is 
found across proto-cluster scales over ~15000 AU, i.e. a 
factor of 5 larger than the linear scales over which the 
chemical segregation toward AFGL2591 VLA 3 has been 
detected (<3000AU). Therefore, the chemical segrega- 
tion toward AFGL2591 VLA 3 is likely related to the 
physical structure and physical processes taking place in 
the molecular envelope around the VLA 3 source. 

One may also consider that the two main peaks seen 
in Type II and III species could be associated with two 
independent sources coincident with the bright 1.3 mm 
continuum source toward AFGL2591 VLA 3 and the 
weaker secondary peak detected toward the south of 
this source (Figure 3). However, species such as SO2, 
OCS or CH3OH show emission peaks which are located 
more than one beam away from these sources. Further- 
more, the global kinematics of the molecular gas toward 
AFGL2591 VLA 3 fit well with Keplerian-like rotation 
when all molecular emission across the envelope is con- 
sidered (see Figures 4 and 5, and Section 3.3), suggesting 
that this emission likely arises from the same physical 
entity rather than being associated with two different 
sources. 

3.3. Kinematics of the molecular gas toward AFGL2591 

VLA 3 

Figure 4 presents the channel maps of some represen- 
tative lines of Type I (H2S), Type H (SO2 and HC3N) 
and Type HI (CH3OH) species. These maps are su- 
perimposed on the 1.3 mm continuum emission of the 
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Fig. 4. — Channel maps of the molecular emission of several lines of H2S (Type I species), SO2 and HC3N (Type II), and CH3OH 
(Type III), observed toward AFGL2591 VLA 3 (solid contours). These maps are superimposed on the 1.3mm continuum emission detected 
toward this source (in colour; dashed contour as in Figure 3). The first contour and step levels of the molecular emission are, respectively, 
0.075 (3a) and 0.075 Jy/beam for H2S (2i,o^2i,i), 0.06 (3a) and 0.06Jy/beam for SO2 (115,7^124,8), 0.06 (3o-) and 0.06Jy/beam for 
HC3N ^;=0 (24^23), and 0.06 (3o-) and 0.06 Jy/beam for CH3OH vt=0 (4+2,2^3+1,2 E). The radial velocities of the channel maps are 
shown at the upper left corner in every panel. The central radial velocity of the AFGL2591 source is vlsr=-^-^^^^~^ (van der Tak et 
al. 1999). Dashed lines show the direction of the PV cuts reported in Figure 5. Symbols are as in Figure 3. The beam size is shown at the 
lower right corner in the first panels from the left. 



AFGL2591 VLA 3 source. From these maps, we find 
that the molecular gas toward VLA 3 shows a velocity 
gradient from the north-east to the south-west of this 
source, from blue- to red-shifted velocities. A similar be- 
havior was reported for HDO and H2^0 by van der Tak 
et al. (2006), and contrasts with that of the outflowing 
gas (Section 3.1). Since the outflow axis lies almost along 
the line-of-sight (van der Tak et al. 1999), the kinematics 
of the molecular gas toward AFGL2591 VLA 3 are con- 
sistent with counter-clockwise rotation (van der Tak et 
al. 2006). We note that the channel maps do not show 
any clear molecular feature arising from the recently re- 
ported source NE (Sanna et al. 2012). 

The presence of a velocity gradient for the molecular 
gas across AFGL2591 VLA 3 is also evident from the 
Position- Velocity (PV) diagrams obtained along the PV 
cuts shown in Figures 3 and 4 (see Figure 5) . The direc- 
tion of the PV cuts (P.A.~20°) has been selected so that 
it is roughly perpendicular to the axis of the AFGL2591 
VLA 3 outflow (P.A.-280°; van der Tak et al. 1999). The 
offsets of the PV diagrams (in arcseconds) are calculated 
from the northeast to the southwest of AFGL2591 VLA 
3. 

From Figures, we find that Type I, II and III species 
show significant differences in their kinematic distribu- 



tion with holes in the PV diagrams of SO2, SO and 
CH3OH. The kinematic structure of all these molecules, 
however, seems to be complementary, as expected if they 
probed different regions within the envelope. As shown 
in Figures, H2S (Type I) traces the regions closer to the 
VLA 3 protostar and is detected across a larger veloc- 
ity range (from -12kms~"^ to Okms""*^) than SO2 (from 
-lOkms"^ to -2kms~^) or CH3OH (from -8kms~^ to 
-3kms~^). In contrast to H2S, SO2 (Type II species) 
peaks at regions further away (at ^0.25"-0.35") from the 
position of the central source, while CH3OII (Type III) 
is found at ^0.5" away from VLA 3. 

Since every species probes a particular region within 
the AFGL2591 VLA 3 hot core, the global gas kinemat- 
ics of the envelope can be obtained by combining the data 
from Type I, Type II and Type III species. For this, we 
have used the molecular line transitions H2S (22 o^2i 1), 
SO2 (115,7^124,8) and CH3OH ^,=0 (4+2,2^3+1,2 E). 
The combined image is generated by normalizing the in- 
dividual hue datacubes to their measured peak intensi- 
ties and, after normalization, the individual images are 
added by considering that each of them accounts for 1/3 
of the final image. The combined image is shown in panel 
h of Figures. We note that if H2S were not included in 
the combined image, the final PV diagram would show 
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Fig. 5. PV diagrams of several lines of H2S, SO2, SO and CH3OH observed toward AFGL2591 VLA 3 (panels from a to g). The 
first step and contour levels are 50 {2a) and 50mJy/beam for H2S (22,0— ^2i^i), 40 (2cr) and 40mJy/beam for SO2 (115,7— )^124^8)5 60 
(3cr) and 60mJy/beam for SO2 (227,15^236, is), 60 (3cr) and 60 mJy /beam 'for SO (56^45), 40 {2a) and 40mJy/beam for CHsOH 
^^=0 (4+2,2^3+1,2 E), 50 {2a) and 25mJy/beam for CH3OH vt=l (61,5^72,6 A"), and 50 {2a) and 25mJy/beam for CH3OH vt=0 
(20+1,19— ^►20+0, 20 E). Offsets (in arcseconds) are calculated along the PV cuts shown in Figures 3 and 4 (P.A.~20°) from the northeast to 
the southwest of AFGL2591 VLA 3. Horizontal and vertical dotted lines indicate the position and radial velocity of the AFGL2591 VLA 
3 source. Panel h) shows the combined map of the H2S (22,0— ^2i,i), SO2 (115,7— )'124,8) and CH3OH vt=0 (4+2,2—^3+1,2 E) lines. The 
first contour and step levels for the combined map are 70 {2a) and 105mJy/beam. Dashed, solid and dotted lines show the PV curves 
calculated assuming Keplerian rotation for a source with central mass 10 M0, 40 and 100 M0, respectively. The envelope size and 
inclination angle assumed for these PV curves are 5400 AU and i=20° , respectively. 



a central hole in its distribution preventing us to char- 
acterize the very inner regions with the most blue- and 
red-shifted velocities toward the AFGL2591 VLA 3 hot 
core. 

In Figures (panel h), we also compare the combined 
map with the expected PV diagrams from a rotating 
envelope that follows a Keplerian law around a central 
source with a mass of 10 M©, 40 and 100 (see 
dashed, solid and dotted lines, respectively). The kine- 
matics of the molecular gas toward the AFGL2591 VLA 
3 envelope can be explained by Keplerian-like rotation 
around a 4OM0-source. The derived envelope size is 
5400 AU, and the inclination angle is i=20°. This an- 
gle is consistent with that proposed by van der Tak et al. 
(1999) and van der Wiel et al. (2011) for this source. The 
derived mass of the central object, 40 M©, is consistent 
with that estimated by Sanna et al. (2012, of ~38M0), 
and corresponds to an 06 star on the ZAMS. Previous 
estimates of the mass of the AFGL2591 VLA 3 source 
were underestimated due to the closer distance consid- 
ered (of '^1 kpc; see e.g. van der Tak et al. 1999; Trinidad 
et al. 2003). 

3.4. UV coverage and missing flux in the VEX images 

The differences in the morphology between Type I, II 
and III species toward AFGL2591 cannot be attributed 



to differences in the UV plane coverage, because all 
molecular lines were observed simultaneously within the 
same track (see Figure 1). 

Alternatively, one may think that the observed mor- 
phological differences are due to missing flux in the VEX 
observations. In Figure 6, we compare the combined 
images of the 1.3 mm continuum, SO2 (Type II) and 
CH3OH emission (Type III species) obtained with the 
VEX + Subcompact data, with those generated with 
the VEX data only. From Figure 6, we find that the 
combined VEX + Subcompact maps provide a very sim- 
ilar spatial distribution to that reported in the VEX im- 
ages. This demonstrates that the observed differences 
in the morphology of Type I, II and III species are not 
produced by missing short spacings, but due to a real 
chemical segregation within the AFGL2591 VLA 3 hot 
core. We note that in Figure 6 we do not include the 
combined image of the H2S (22,o^2i,i) transition, be- 
cause this line was not covered within the Subcompact 
observations (see Section 2). 

4. EXCITATION TEMPERATURE, COLUMN DENSITIES 
AND MOLECULAR ABUNDANCES TOWARD AFGL2591 

From the integrated intensities of the SO2, CH3OH 
and HC3N lines measured toward single beams across 
the AFGL2591 VLA 3 hot core, we can estimate the 
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Fig. 6. — Comparison of the combined (VEX + Subcompact) in- 
tegrated intensity images of the 1.3 mm continuum emission, and of 
the SO2 (115,7^124,8) and CH3OH vt=0 (4+2,2^3+1,2 E) hues, 
with those obtained by using the VEX data only. The velocity 
range considered to generate the SO2 and CH3OH images goes 
from -8.7 to -2.1kms~^. For the 1.3mm images (color and thin 
contours in left panels), the contour levels are 6.0 (4<t), 7.5, 10.5, 
25.5, 40.5 and 55.5mJy/beam for the VEX + Subcompact data, 
and of 3.2 (4a), 4.0, 5.6, 17.6, 29.6, 41.6 and 53.6 mjy/beam for the 
VEX map. The SO2 and CH3OH emission maps (thick contours 
in central and right panels) appear superimposed on the 1.3 mm 
continuum images (color scale and dashed contour). For the SO2 
and CH3OH line data, the first contour and step levels are, respec- 
tively, 0.5 (5<t) and 0.3 Jy/beam km s~^ for the combined VEX + 
Subcompact images, and 0.35 {5a) and 0.21Jy/beam kms~^ for 
the VEX maps. Symbols are as in Figures 3 and 4. Beams are 
shown at the lower right corner in every panel. 



excitation temperature of the gas across the envelope, 
Tex, by means of the population diagram method that 
assumes optically thin emission and LTE (Goldsmith & 
Langer 1999). We select the Continuum Peak position, 
and the emission peaks Tl and T2 because they show 
the brightest molecular emission for Type I, II and III 
species, respectively (see Figure 3). In Table 2, we report 
the observed parameters of the single beam lines mea- 
sured toward the Continuum Peak for H2S and ^^CS, 
and toward emission peaks Tl for SO, OCS, HC3N and 
SO2, and T2 for CH3OH (see Figure 3 and this Table 
for the exact coordinates of these positions). The de- 
rived central radial velocities, vlsr^ ^ind linewidths, Af , 
typically range from -5.5 to -6.5kms~^, and from ^^2.5 
to 4kms~^, respectively, although significant discrepan- 
cies in vlsr and Av are found for some of the species 
within the same type (see e.g. SO for Type II or the 
low-excitation lines CH3OH Vt=0 (5+i 4^4+2 2 E) and 
CH3OH vt=0 (3_2,2^4_i,4 E) for Type III; ' Table 2). 
These discrepancies could be due to large optical depths 
associated with these low-excitation lines {Eu<100 K; see 
Table 1). 

In Figure?, we show the population diagrams of SO2 
and HC3N derived toward position Tl, and the popu- 
lation diagram of CH3OH obtained toward position T2. 
The derived Tex are 185±11K for SO2, 150±20K for 
HC3N, and 124±12K for CH3OH, suggesting a temper- 
ature gradient within the envelope. The derived column 
densities of SO2, HC3N and CH3OH are 2xl0i^cm-2, 
8xl0^^cm~^ and 4xl0^^cm~^, respectively. In the case 
of CH3OH, was derived from lines with ^^i>300K 
because its lower-excitation transitions are likely opti- 



cally thick. As shown by Girart et al. (2002), large op- 
tical depths indeed lead to the underestimation of the 
measured column densities of the low-excitation transi- 
tions of CH3OH, clearly deviating from the optically thin 
solution provided by the population diagram method. 

For H2S, -^^CS, OCS, and SO, the molecular column 
densities are estimated from single transitions by consid- 
ering optically thin emission and LTE conditions. The 
assumed excitation temperature for all these lines is 
Tea; =185 K, the excitation temperature found for SO2 to- 
ward Tl. For ^^CS, we have also considered an isotopic 
ratio I2c/i3c=60 (Wilson & Rood 1994). The derived 
column densities are 7x10^^ cm~^ for H2S, 4x10^^ cm~^ 
for i^CS, 2xl0i^cm-2 for OCS, and SxlO^^cm-^ for 
SO. 

Since H2S and ^^CS are mainly detected toward the 
Continuum Peak, it is possible that their excitation tem- 
perature is higher than assumed here, as suggested by our 
modelling of the chemical segregation toward AFGL2591 
VLA 3 (see Section 5). In case Te^ were ~1000K (as 
measured by Carr et al. 1995, from ro-vibrational ab- 
sorption fines of HCN), tfie column densities derived for 
H2S and ^^CS would be factors ~5-8 fiigfier tfian tfiose 
reported above. 

To estimate the molecular abundances across 
the AFGL2591 VLA 3 envelope, we used the 
1.3 mm continuum fiux toward tfie Continuum Peak 
(68 mJy/beam), toward Tl (19 mJy/beam), and toward 
T2 (8mJy/beam), to derive tfie H2 column densi- 
ties at tfiese positions. From Equation 2 in Enocfi 
et al. (2006), and by considering dust opacities of 
O.Olcm^g"^ (Ossenkopf & Henning 1994), we estimate 
H2 column densities of 2xl0^^cm~^, 4xl0^^cm~^ and 
3xl0^^cm~^ toward tfie Continuum Peak, Tl and 
T2, respectively. We fiave assumed dust temperatures 
of -200 K toward tfie Continuum Peak, of -200 K 
toward Tl, and of — 120K toward T2, following tfie 
gas excitation temperatures derived from SO2 and 
CH3OH (see above). However, we note tfiat tfiese 
H2 column densities are subject to uncertainties since 
tfie dust temperature cannot be constrained from our 
observations. If tfie temperature of tfie dust were fiigfier 
tfian 200 K toward tfie Continuum Peak as suggested by 
tfie modeling of van der Tak et al. (1999, -500-800 K for 
angular scales <0.2"), tfie derived H2 column densities 
would be decreased to — 3-5x10^^ cm~^. 

Tfie derived molecular abundances are presented in 
Columns 2 and 4 of Table 3. In tfiis Table, we also in- 
clude tfie HCN abundances measured toward AFGL2591 
and reported by Doty et al. (2002). HCN is detected 
in absorption in tfie IR toward tfiis source (Carr et al. 
1995); and since absorption is restricted to tfie direction 
of tfie line-of-sigfit toward tfie source, tfie bulk of tfie 
HCN gas is likely associated witfi tfie innermost regions 
of tfie AFGL2591 VLA 3 envelope. We do not include 
tfie abundance from otfier molecular species reported in 
Doty et al. (2002), because tfiey are mainly obtained from 
emission measurements toward tfiis source, wfiicfi likely 
probe gas at all scales tfirougfiout tfie envelope. 

From Table 3, we find tfiat tfie derived abundances de- 
crease witfi increasing radius from tfie position of tfie 
Continuum peak to positions Tl and T2 by factors of 
>5-10 for H2S and CS, and by two orders of magnitude 
for HCN. For HC3N, SO2, OCS and SO, fiowever, tfieir 
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TABLE 2 



Molecule 


Transition 




Continuum 


Peak^ 






'^LSR (kms ^) 


/\v (kms~-^) 


Peak Intensity (Jy/beam) 


H2S 


22,o^2i^i 


-5.4(0.5) 


7.3(1.2) 


0.34(0.03) 




'5^4 ' 


-5.44(0.16) 


3.7(0.5) 


0.18(0.02) 








Position Tl"^ 






n ^T\— 


Ai; (kms~-^) 


Peak Intensity (Jy/beam) 


Qn 


06—^45 


A 71 ( C\ C\Ci\ 

-4. ( i(^U.Uo j 


6.42(0.13) 


0.60(0.02) 


U^D 




-D.z / (u.uy j 


3.6(0.2) 


0.352(0.016) 




1 Q v1 Q 


-D.zH^U. 10 j 


4.5(0.4) 


0.21(0.02) 




n o/i V OQ 
V=U z4 — >-Zo 


-D.o(U.i ) 


3.7(0.2) 


0.28(0.02) 




1 O/l vOQ 

i;7 — ie ^4 — ?-ZO 


1 o/'n /I \ 
- / .Z(^U.4 j 


1.3(0.0) 


0.10(0.03) 




lyy — ±+ ^"-t r^Kj 


-7.7(0.4) 


1.7(0.5) 


0.06(0.02) 




VQ=le 24^23 






<0.06^ 




i;6=l/ 24^23 






<0.05 


SO2 


227,15— )-236, 18 


-5.81(0.12) 


3.1(0.3) 


0.193(0.017) 




222,20^22i,2i 






<0.05 




163,13^162,14 






<0.05 




115,7^124,8 


-5.69(0.07) 


3.41(0.18) 


0.356(0.018) 




132,12^13l,l3 






<0.06 




37io, 28^389, 29 






<0.05 




143,11— >142, 12 






<0.05 



Position T2' 





VLSR (kms ^) 


(kms ^) 


Peak Intensity (Jy/beam) 


CH3OH ^£=0 5+1,4^4+2,2 E 


-5.3(0.5) 


4.1(1.0) 


0.07(0.03) 


-yt^l 61,5^72,6 A- 


-6.41(0.13) 


2.4(0.3) 


0.174(0.018) 


vt=^ 20+1,19^20+0,20 E 


-5.95(0.15) 


2.1(0.4) 


0.18(0.02) 


Vt=0 4+2,2^3 + 1,2 E 


-5.93(0.07) 


3.29(0.15) 


0.39(0.02) 


vt=0 8+0,8^7+1,6 E 


-5.7(0.3) 


4.9(0.7) 


0.087(0.017) 


Vt=0 15+4,11^16 + 3,13 E 


-6.0(0.3) 


4.4(0.8) 


0.085(0.018) 


vt=0 8-1,8^7+0,7 E 


-5.74(0.08) 


3.14(0.18) 


0.290(0.018) 


Vt=0 195,15-^204,16 A+ 






<0.05 


vt=0 195,14^204,17 A- 






<0.04 


Vt=0 3-2,2^4-1,4 E 


-7.13(0.15) 


2.5(0.4) 


0.124(0.016) 


vt=0 102,9^93,6 A- 


-6.63(0.17) 


3.8(0.4) 


0.19(0.02) 


vt=0 102,8^93,7 


-6.3(0.2) 


4.1(0.5) 


0.16(0.02) 



^The coordinates of the Continuum Peak are a( J2000)=20^29"^24.889^ (5(J2000)=40° 11^19.51" . 
^The coordinates for the Tl position are a( J2000)=20^29'^24.870^ (5( J2000)=40° 11'19.73" . 
^The coordinates for position T2 are a( J2000)=20^29"^24.884^ (5(J2000)=40° 11^20.00" . 
^Upper limits refer to the 3a level in the extracted single beam line spectra. 



derived abundances are factors >6 larger toward posi- 
tion Tl than toward the Continuum Peak. Fohowing a 
similar trend, the derived abundance of CH3OH is a fac- 
tor of >100 larger toward position T2 than toward the 
Continuum Peak. 

We note that, if the assumed excitation temperature 
of the gas were ^800 K toward the continuum peak, 
the measured molecular abundances toward this position 
(see Column 2 in Table 3) would be increased by factors 
of ~2-40. This implies that the increasing abundance 
trend measured for Type II species (i.e. HC3N, SO2, 
OCS and SO) across AFGL2591 VLA 3 would become 
non-existent. However, in the case of CH3OH (Type 
III species), the derived upper limit to its abundance 
(<2xl0-^) would stih be one order of magnitude lower 
than the CH3OH abundance measured toward position 
T2 (Table 3); and for Type I species (H2S and CS), their 
decreasing abundance trend with larger radii would be 
even stronger than reported in Table 3 (the H2S and CS 
abundances measured toward the Continuum Peak would 
be factors of ~ 100-500 higher than toward T2). 

We finally stress that the scope of this paper is to show 
the abundance trends found across the AFGL2591 VLA 
3 envelope, and not to obtain the global molecular abun- 
dances measured toward this object for which observa- 
tions with lower angular resolution would be required. 



5. CHEMICAL MODELING OF THE AFGL2591 HOT CORE 

To qualitatively explain the chemical segregation ob- 
served toward AFGL2591 VLA 3, we have used the two- 
phase UCL CHEM gas-grain chemical code (Viti et al. 
2004). In Phase I, this code simulates the free-fall grav- 
itational collapse of a cloud where the mantles of dust 
grains form via hydrogenation reactions. The initial 
atomic abundances are solar (Sofia & Meyer 2001). The 
formation of the initial core occurs from a diffuse medium 
with H2 densities of ~100cm~^ at a kinetic tempera- 
ture of 20 K. The core collapse is stopped when the fi- 
nal density is reached in the model, and it corresponds 
to time-scales of ~5xl0^yrs. In Phase 2, we calculate 
the time-dependent evolution of the chemistry of gas and 
dust once stellar activity is present. In our case, we as- 
sume that the mantles are instantaneously evaporated 
after the turning on of the protostar. We also assume 
that H2S is the main reservoir of atomic sulfur on dust 
grains. 

The detection of ro- vibrational absorption of HCN to- 
ward AFGL2591 VLA 3 (Carr et al. 1995; Lahuis & van 
Dishoeck 1997), suggests the presence of an inner and 
hotter core with relatively low extinction (^^^30"^), em- 
bedded within a more extended and cooler envelope (van 
der Tak et al. 1999). The AFGL2591 VLA 3 hot core 
has therefore been modeled taking into account these 
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Fig. 7. — Population diagrams for SO2, HC3N and CH3OH derived from the line spectra extracted from positions Tl and T2 toward 
AFGL2591 VLA 3 (see Table 2). The derived Tex are shown in the upper right corner of every panel. Error bars correspond to la 
uncertainties and vertical arrows indicate upper limits. 



TABLE 3 


Abundances 


Observed 


ModeP 


Observed 


ModeP 




(Continuum Peak) 


(Region A) 


(Positions T1/T2) 


(Region B) 


H2S 


5xl0-« 


5x10-9 


<4xlO-9b 


1x10-10 


CS 


2x10-^ 


1x10-6 


<4xl0-9b 


1x10-^ 


HON 


1x10-6^ 


4x10-^ 


1x10-8^ 


3x10-^ 


HC3N 


7x10-9 


4x10-11 


2x10-^ 


3x10-^ 


OCS 


9x10-9 


1x10-11 


6x10-^ 


2x10-9 


SO2 


1x10-^ 


3x10-11 


6x10-^ 


1x10-6 


SO 


1x10-^ 


2x10-9 


6x10-^ 


1x10-^ 


CH3OH 


<2xl0-^^ 


<lxlO-i^ 


2x10-6^ 


5x10-^ 



^Abundances predicted by the model for a dynamical age of ~2xl0^yr (Doty et al. 2002). 
^Upper limits derived from the la level in the H2S and i^CS spectra toward position T2. 
^From Doty et al. (2002). 

•^Upper limits derived from the la level in the CH3OH spectrum toward the position of the Continuum Peak. 
^CH30H abundance derived from the population diagram shown in Figure 7 for position T2. 



two physical regimes (regions A and B; Table 4). Re- 
gion A corresponds to the inner core where the emission 
from H2S and ^^CS peaks (i.e. toward the Continuum 
Peak with angular scales of '-^0.2" or r=600 AU at a dis- 
tance of 3 kpc). Region B corresponds to an intermediate 
position between Tl and T2 in our observations where 
HC3N, OCS, SO2, SO and CH3OH show their maximum 
emission (i.e. at ~0.4", or r=1100 AU, from the central 
source in our images; see Figure 3). Therefore, for sim- 
plicity, our two region model does not discriminate be- 
tween positions Tl and T2 within the AFGL2591 VLA 3 
envelope. We also consider the existence of an inner cav- 
ity within the AFGL2591 VLA 3 hot core, as proposed 
by Preibisch et al. (2003) and de Wit et al. (2009). These 
authors derived a radius of 40 AU for this cavity assum- 
ing a distance of ~1 kpc for the AFGL2591 VLA 3 source 
(van der Tak et al. 1999). Since the actual distance to 
this source is 3 kpc (Rygl et al. 2012), the radius of the 
cavity considered in our model is 120 AU, i.e. a factor of 
3 larger than that derived by Preibisch et al. (2003) or 
de Wit et al. (2009). This cavity remains unresolved in 
our VEX images. The wall of the cavity at r=120 AU is 
the position at which the visual extinction in our model 
is Ay=0. 

The H2 densities at the end of Phase I for regions A 
and B are estimated by using the H2 density distribu- 
tion calculated by van der Tak et al. (2000a), after cor- 
recting it from the distance of 3 kpc (see discussion in 
Section 4.2 of van der Tak et al. 1999, on how the H2 
density distribution for AFGL2591 VLA 3 is affected by 
a larger distance to the source). The averaged H2 densi- 
ties for regions A and B are, respectively, 4.6xl0^cm~^ 
and 1.8xl0^cm~^. These densities are kept constant in 



Phase II of our model. 

In order to compare the H2 column densities expected 
from these H2 density values with those measured to- 
ward the Continuum Peak and positions T1/T2, one 
needs to consider the total H2 column density across the 
AFGL2591 VLA 3 envelope. For instance, to compare 
the results from region A with those observed toward the 
Continuum Peak, the contribution from the inner core 
(-5xl0^cm-^x960 AU, i.e. 2x480 AU, the shell size for 
region A; see below) must be added to that from the outer 
envelope [-2 x 10^ cm-^ x (5400-960) AU, where 5400 AU 
is the envelope size derived in Section 3.3]. This gives a 
total H2 column density of ~2xl0^^cm~^. For region B 
(or positions T1/T2), the total H2 column density can 
be estimated simply by doing ~2xl0^ cm~^x5400 AU, 
which gives ~lxl0^^cm~^. For region B, the expected 
and observed H2 column densities are within a fac- 
tor of ^4. However, for region A (or the Continuum 
Peak), the low value of the expected H2 column density 
(2xl0^^cm~^ vs. 2xl0^^cm~^; see Section 4) suggests 
that the actual temperature of the dust toward the Con- 
tinuum Peak is likely higher than 200 K (Section 4). 

The gas and dust temperatures assumed for regions A 
and B are obtained from the theoretical modelling of the 
core by Stauber et al. (2004) for region A (-1000 K), 
and by van der Tak et al. (1999) for region B (-200 K). 
We stress that the dust and gas temperature profiles do 
not need to be corrected by the distance of 3 kpc because 
they are constant in terms of projected distance (i.e. in 
arcseconds; see Section 4.2 of van der Tak et al. 1999). 

To derive the FUV luminosity emitted by the 
AFGL2591 VLA 3 source, we have used the FUV lu- 
minosity profile derived by Bruderer et al. (2009a) as a 
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Fig. 8. Abundances of H2S, CS, SO2, SO, HON, HC3N, OCS 
and CH3OH predicted for Phase 2 by the two-point chemical model 
of the AFGL2591 VLA 3 envelope (Sections). The physical prop- 
erties assumed for regions A and B are reported in Table 4. Ver- 
tical dotted lines show the dynamical age of the AFGL2591 enve- 
lope (~2xl0^yr) derived by Doty et al. (2002) and Stauber et al. 
(2005). 



TABLE 4 





Region A 


Region B 


r 


600 AU 


1100 AU 


n(H2)^ 


4.6x10^ cm-3 


1.8x10^ cm-3 


Ti ■ ^ 


1000 K 


200 K 


Av^ 


20.5^ 


28.5^ 



^From Boonman et al. (2001). 

'^Consistent with HCN line absorption data (Carr et al. 1995; van 
der Tak et al. 1999) and with our SO2 data (Section 4). 
•^See text for the derivation of Ay. 



function of effective temperature, Te// (see Figure A. 1 in 
their work). For an 06-type star on the ZAMS with an 
effective temperature of Te//~4xlO^K, the luminosity 
in the FUV band is ~2 x 10^^ erg s~ ^ per . Considering 
that the AFGL2591 VLA 3 source has a total bolometric 
luminosity of L5o^~2x 10^ (Section 3 and Sanna et al. 
2012), we obtain that the total FUV luminosity emit- 
ted by this source is 4x 10^^ ergs~^. This luminosity is 
a factor of 10 higher than that used by Bruderer et al. 
(2009a), because the distance assumed in their calcula- 
tions was 1 kpc instead of 3 kpc. Since our model consid- 
ers a cavity within the AFGL2591 VLA 3 envelope, the 
FUV photon flux illuminating the wall of this cavity (at 
120 AU) is -6xl0^hab (1 hab=1.6xl0-^ergs-^ cm-^; 
Habing 1968), as derived from geometrical dilution. We 
thus assume an incident FUV flux of '^GxlO^hab at the 
position with A^=0 in our model. 

At r>120AU, the extinction is calculated as 
A^ = [sizexn(H2)]/[1.6xl0^^cm~^]. For region A (sheh 
thickness 5zze=600 AU-120 AU=480 AU), the extinction 
is ~20.5^. For region B, however, the individual ex- 
tinction from regions A and B should be added since 
the FUV field from the AFGL2591 star wih be atten- 
uated as it propagates outwards through the envelope. 
The individual extinction at r=1100 AU (shell thickness 
5ize-1100 AU-600 AU=500 AU) is -8^, leading to a to- 
tal extinction of -28.5^ (20.5^+8^) for region B. 

In Figure 8, we present the abundances of H2S, CS, 
HCN, SO, SO2, HC3N, OCS and CH3OH, predicted 
for Phase 2 by our two-point model of the AFGL2591 



VLA 3 envelope. Our results show that the chemi- 
cal segregation observed toward this source can be ex- 
plained by the combination of molecular UV photo- 
dissociation and a high-temperature gas-phase chemistry 
in the inner core of AFGL2591 VLA 3. In region A, 
H2S and CII3OH are initially destroyed by FUV pho- 
tons (extincted FUV flux of -7.5 hab for A^-20.5^) 
after the release of the grain mantles into the gas- 
phase (t<100yr). However, the high temperatures of 
the gas (T/ein'^lOOO K) allow the re-formation of H2S 
via the endothermic reaction H2+HS^H2S + H (acti- 
vation barrier of 8050 K), making H2S abundant in the 
inner core (5xl0~^ for t~2xl0'^yr; Figure 8). We note 
that these high temperatures are strictly required for re- 
gion A because the H2S abundance would dramatically 
drop otherwise, in contrast with what is observed. Like 
H2S, HCN and CS survive the FUV photo-dissociation 
thanks to the gas-phase reactions H2 + CN^HCN + H 
and C2 + S^CS-hC. CH3OH is completely destroyed 
since it cannot be re-formed in the gas phase (e.g. Garrod 
et al. 2008). 

In region B, molecular photo-dissociation is less severe 
(extincted FUV flux of only -0.003 hab for A^-28.5^); 
and since the gas temperature is only 200 K, the ac- 
tivation barrier for the endothermic reaction to form 
H2S cannot be overcome. As a consequence, other S- 
bearing products such as SO, OCS and SO2 are enhanced 
(Figure 8). This explains the detection of H2S only to- 
ward the inner core of AFGL2591 and the double-peaked 
structure observed for SO, SO2 and OCS (Section 3). 
We note that this S-bearing chemical structure was also 
predicted by more complex chemical modeling of the 
AFGL2591 VLA 3 source, which includes UV and X-ray 
photo-chemistry (Bruderer et al. 2009b). 

In the case of CH3OH and HC3N, these molecules are 
found to survive in region B thanks to the higher ex- 
tinction. For HCN, the formation of this molecule is less 
efficient at lower temperatures, explaining the decrease 
of the HCN abundance (by a factor of >100) in region 
B (Figure 8). This decrease is consistent with that ob- 
served toward the AFGL2591 VLA 3 envelope (Table 3; 
van der Tak et al. 1999; Boonman et al. 2001). 

The agreement between the predicted and measured 
abundances for the rest of the molecular species is within 
factors —2-50, except for HC3N, OCS and SO2 for region 
A where the discrepancies are found to be of several or- 
ders of magnitude. This could be due to the simplistic 
physical structure assumed for the AFGL2591 VLA 3 
envelope. Continuous variations of the density and tem- 
perature of the gas and dust with increasing radius, as 
derived by comprehensive radiative transfer and chemi- 
cal models (see e.g. van der Tak et al. 1999; Doty et al. 
2002; Bruderer et al. 2009b), could significantly change 
the predicted abundances of these molecules. Indeed, 
the molecular abundances are very sensitive to the ex- 
tinction assumed for regions A and B {Ay in Table 4), 
which depends on the H2 density distribution of the en- 
velope and the distance to the central protostar. If the 
assumed extinction is decreased by 1^ for regions A and 
B, the general trend for the molecular abundance is to 
decrease by factors of —3-7, except for HC3N, H2S and 
OCS in region B where their abundances are decreased 
by factors >20. 

Alternatively, the large discrepancies in the measured 
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and predicted abundances of HC3N, OCS and SO2 for re- 
gion A could be due to the fact that the beam size of our 
VEX observations is larger (0.3" xO.5") than the angular 
scale considered for region A (of only 0.2"). The mea- 
sured abundances for HC3N, OCS and SO2 could there- 
fore have received some contribution from emission out- 
side the inner 0.2" of the envelope, artificially increasing 
the abundances for these molecules toward this region. 

We note that the maximum extinction allowed for re- 
gion B is 28.5"^. If an extinction of Ay>29'^ is assumed 
for this region, the abundance of H2S increases to ~10~^, 
which is several orders of magnitude higher than found 
from our observations (the upper limit to the H2S abun- 
dance measured toward position T2 is of <4xl0~^; see 
Tables). 

Finally, our results do not depend on our selection of 
H2S as the main carrier of S on the grain mantles. We 
have run two models with i) 100% of S locked onto OCS; 
and h) 50% of S locked on H2S and 50% just as S. For 
region A, once OCS is released from the mantles, it is 
photo-dissociated leading to large abundances of S. S 
subsequently reacts to form HS, which is rapidly con- 
verted into H2S. This also applies to the case with 50% 
S and 50% H2S, where all S is transformed into H2S. For 
region B, the only difference is found for the 100% OCS 
model, where OCS is the most abundant molecule after 
the release of the mantles. OCS is then efficiently con- 
verted into SO and SO2, showing a similar behavior to 
that found for our previous 100% H2S model. 

6. CONCLUSIONS 

We report the detection of a clear chemical segrega- 
tion at linear scales <3000AU within the massive hot 



core around the high-mass protostar AFGL2591 VLA3 
with the SMA. The molecular gas is distributed in three 
concentric shells where: i) H2S and ^^CS (Type I species) 
peak at the AFGL2591 VLA 3 continuum source; ii) 
HC3N, OCS, SO2 and SO (Type II species) show a 
double-peaked morphology circumventing the dust con- 
tinuum peak; and iii) CH3OH (Type III species) is dis- 
tributed in a ring-like structure that surrounds the dust 
continuum emission. The global kinematics of the gas, 
probed by the different molecules found across the enve- 
lope, are consistent with Keplerian-like rotation around 
a central source of 40 M©. This mass is consistent 
with previous findings (Sanna et al. 2012). In addi- 
tion, the molecular gas shows a temperature gradient 
for increasing distance to the central source, as revealed 
by the higher Tex derived for SO2 (of 185=bllK) than 
for CH3OH (of 124±12K). By comparing our SMA im- 
ages with gas and dust chemical modeling, we find that 
the chemical segregation toward the AFGL2591 VLA 3 
hot core is explained by the combination of molecular 
UV photo-dissociation and a high-temperature gas-phase 
chemistry within the low extinction innermost region in 
the AFGL2591 VLA 3 massive hot core. 
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